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Abstract 

Effective-field-theory methods are used to separate the free energy for a nonabelian 
gauge theory at high temperature T into the contributions from the momentum scales 
T, gT, and g 2 T, where g is the coupling constant at the scale 2ttT. The effects of the 
scale T enter through the coefficients in the effective lagrangian for the 3-dimensional 
effective theory obtained by dimensional reduction. These coefficients can be calculated 
as power series in g 2 . The contribution to the free energy from the scale gT can be 
calculated using perturbative methods in the effective theory. It can be expressed as 
an expansion in g starting at order g 3 . The contribution from the scale g 2 T must 
be calculated using nonperturbative methods, but nevertheless it can be expanded in 
powers of g beginning at order g e . We calculate the free energy explicitly to order g 5 . 
We also outline the calculations necessary to obtain the free energy to order g 6 . 



1 Introduction 



One of the most dramatic predictions of quantum chromodynamics (QCD) is that when 
hadronic matter is raised to a sufficiently high temperature or density, it will undergo a 
phase transition to a quark-gluon plasma. One of the major thrusts of nuclear physics in the 
next decade will be the effort to study the quark-gluon plasma through relativistic heavy-ion 
collisions. For this effort to be successful, it will be important to understand the properties 
of the plasma as accurately as possible. The two major theoretical tools that have been used 
to study the quark-gluon plasma are lattice gauge theory and perturbative QCD. Lattice 
gauge theory has the advantage that it is a nonperturbative method and applies equally well 
to the quark-gluon phase and to the hadron phase. It is an effective method for calculating 
the static equilibrium properties of hadronic matter with baryon density. Unfortunately, 
the Monte Carlo methods used in lattice gauge theory can not be easily applied to problems 
involving dynamical properties or to hadronic matter that is away from thermal equilibrium 
or has nonzero baryon density. These are severe restrictions, because a quark-gluon plasma 
that is produced in heavy-ion collisions will not be at thermal equilibrium and it may have 
nonzero baryon density. Furthermore, many of the most promising signatures for a quark- 
gluon plasma involve dynamical properties. 

Perturbative QCD can help fill this gap, at least for the quark-gluon phase of hadronic 
matter. This method can certainly be applied to the static equilibrium properties of a quark- 
gluon plasma at baryon density, but there are no apparent obstacles to also applying it to 
dynamical problems, or to non-equilibrium situations, or to a plasma with nonzero baryon 
density. Thus it is a powerful tool for studying various aspects of the quark-gluon plasma 
that might be probed through heavy-ion collisions. However there are potential difficulties 
in applying perturbative QCD to the quark-gluon plasma. The method is based on treating 
the coupling constant g as a small parameter, but g{ji) is a parameter that varies rather 
dramatically with the momentum scale /i. In order to apply perturbative QCD, it is necessary 
that g be small at the scale of the typical momentum of a particle in the plasma, which is of 
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order T or perhaps 2irT. While this is necesssary, it may not be sufficient. At sufficiently high 
order in perturbation theory, any observable becomes sensitive to low momentum gluons that 
interact with a large coupling strength g. In order to rigorously apply perturbative QCD, it 
is essential to be able to unravel the various momentum scales that play an important role 
in a problem. If low- momentum contributions are important, they must be treated using 
nonperturbative methods. 

For a quark-gluon plasma at high temperature, there is a hierarchy of 3 momentum scales 
that play an important role in static properties. First, there is the scale T of the typical 
momentum of a particle in the plasma. Next, there is the scale gT associated with the 
screening of color-electric forces by the plasma. Finally, there is the scale g 2 T associated 
with color-magnetic screening. Only recently has a method been developed that can sys- 
tematically unravel the contributions from these various momentum scales. The method 
is based on the construction of effective field theories that reproduce static observables at 
successively longer distance scales. This effective-field-theory approach is based on an old 
idea called "dimensional reduction" [[I], According to this idea, the static properties of a 
(3+l)-dimensional field theory at high temperature can be expressed in terms of an effective 
field theory in 3 space dimensions. Dimensional reduction has long been used to provide 
insight into the qualitative behavior of field theories at high temperature H], 0, BL M. The 
effective-field-theory approach makes this idea into a practical tool for quantitative calcula- 
tions. In Ref. ||, we developed the effective-field-theory approach to dimensional reduction 
and applied it to a scalar field with a <f) A interaction. We demonstrated the power of this 
method by using it to carry out several perturbative calculations beyond the frontiers set 
by previous work. A similar approach was developed independently by Farakos, Kajantie, 
Rummukainen, and Shaposhnikov B, who applied it to the important problem of the elec- 
troweak phase transition. This method has also been applied to QCD H, and used to resolve 
a longstanding problem involving the breakdown of the perturbation expansion for the free 
energy [g]. These ideas have also been used to determine the asymptotic behavior of the 
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correlator of Polyakov loop operators 0] and to provide a rigorous nonperturbative definition 



of the Debye screening mass in nonabelian gauge theories [1C[. 

Once we have understood how to resolve the contributions of the various momentum 
scales in thermal QCD, asymptotic freedom guarantees us that perturbation theory will be 
under control in the high temperature limit. At sufficiently high temperature, the running 
coupling constant will be small enough that calculations to leading order in g will be accurate. 
However, in most practical applications, such as those encountered in heavy ion collisions, the 
temperature is not asymptotically large, and we must worry about higher order corrections. 
The accuracy of the perturbation expansion can only be assessed by carrying out explicit 
perturbative calculations beyond leading order. One of the obstacles to progress in high 
temperature field theory has been that the technology for perturbative calculations was 
not well developed. Only very recently have there been any calculations to a high enough 
order that the running of the coupling constant comes into play. The simplest physical 
observable that can be calculated in perturbation theory is the free energy, which determines 
all the static thermodynamic properties of the system. The running of the coupling constant 
first enters at order g A . The free energy for gauge theories at temperature but large 
chemical potential was calculated to order g 4 long ago [O. The first such calculation at 
high temperature was the free energy of a scalar field theory with a <j) 4 interaction, which 
was calculated to order g A by Frenkel, Saa, and Taylor in 1992 [fT2"| . (A technical error was 
later corrected by Arnold and Zhai JT3J.) The analogous calculations for gauge theories were 
carried out in 1994. The free energy for QED was calculated to order e 4 by Coriano and 
Parwani [13] and the free energy for a nonabelian gauge theory was calculated to order g A 
by Arnold and Zhai [I3|. The calculation of Arnold and Zhai was completely analytic, and 
thus represent a particularly significant leap in calculational technology. The calculational 
frontier has since been extended to fifth order in the coupling constant by Parwani and 
Singh |15| and by Braaten and Nieto || for the 4 field theory, by Parwani |l(| and by 
Andersen [17j for QED, and by Kastening and Zhai [JlS ] for nonabelian gauge theories. In 
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this paper, we present an independent calculation of the free energy for a nonabelian gauge 
theory to order g 5 |19]|, verifing the result of Kastening and Zhai. In our calculation, we use 
effective-field-theory methods to simplify the calculation and to resolve the contributions to 
the free energy from the momentum scales T and gT. We also outline the calculations that 
are required to obtain the free energy to order g 6 . 

In Section 2, we describe how effective field theories can be used to resolve the contribu- 
tions to the free energy from the momentum scales T, gT, and g 2 T. In Section 3, we calculate 
the coefficients in the lagrangian for the effective field theory obtained by dimensional re- 
duction. In Section 4, we use the effective field theory to calculate the free energy for QCD 
to order g 5 . In Section 5, we outline the calculations that would be necessary to improve 
the accuracy to order g e . In Section 6, we discuss the implications of our calculation for 
convergence of the perturbation expansion for the free energy. We present some conclusions 
in Section 7. In two appendices, we tabulate the analytic expressions for all the sums and 
integrals that arise in our calculation. 



2 Separation of Scales in the Free Energy 

The free energy for QCD at high temperature T includes contributions ^from the momentum 
scales T, gT, and g 2 T. In this section, we explain how the contributions from these three 
momentum scales can be unraveled by using effective field theory methods. 

The static equilibrium properties for hot QCD are given by the free energy density F , 
which is proportional to the logarithm of the partition function: 

F = -^logZ QCD , (1) 

where V is the volume of space. In the imaginary-time formalism for thermal QCD, the 
partition function is given by a functional integral over quark and gluon fields on a 4- 
dimensional Euclidean space. The Euclidean time r is periodic with period (3 = 1/T. The 
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partition function is 

Zqcb = J VA^T)Vq(^T)Vq^r)^ V (^- dr J C QCl )j . (2) 

The gluon fields are periodic functions of r while the quark and antiquark fields are antiperi- 
odic. The lagrangian is 

£qcd = \G%G% + q llx D,q , (3) 

where G a = d^A^ — d u A® + gf abc A b ^A1 is the field strength and g is the gauge coupling 
constant. All the quark fields have been assembled into the multi-component spinor q, and 
the gauge- covariant derivative acting on this spinor is = d^ + igA a ^T a . The relevant quark 
flavors are all assumed to be massless. 

In order to make our calculations as general as possible, we will express them in terms 
of the group-theory factors Ca, Cf, and Tp defined by 

jabcjabd = C A 5 Cd , (4) 

(T a T% = C F S l3 , (5) 
tr (T a T b ) = T F 5 ab . (6) 

For an SU(N C ) gauge theory with nf quarks in the fundamental representation, these factors 
are Ca = N c , Cp = (N% — 1)/(2N C ), and T F = rif/2. The dimensions of the adjoint 
representation and the fermion representation are = N% — 1 and dp = N c rif, respectively. 

The free energy for QCD can also be calculated using an effective field theory in 3 space 
dimensions called electrostatic QCD (EQCD). This effective theory is constructed so that it 
reproduces static gauge- invariant correlators of QCD at distances of order l/(gT) or larger. 
It contains an electrostatic gauge field Ag(x) and a magnetostatic gauge field A"(x) that 
can be identified, up to field redefinitions, with the zero-frequency modes of the gluon field 
y4^(x, r) for thermal QCD in a static gauge 0. The free energy for thermal QCD can be 
written 

F = T (f E (A B ) - , (7) 
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where ^eqcd is the partition function for EQCD: 

^eqcd = PA (x)PA(x) exp (- J d 3 x £ E qcd) • (8) 

The functional integral requires an ultraviolet cutoff A#. The lagrangian for EQCD is 

£eqcd = \G%G% + ^(DiA ) a (DiA ) a + \m\AlA% + \\ E {A a Q Alf + 5£ EQC d, (9) 

where = diA - — djA\ + gEf abc A h i A C j is the magnetostatic field strength with coupling 
constant g E . If the fields A§ and Ai are assigned the scaling dimension 1/2, then the operators 
shown explicitly in (P) have dimensions 3, 3, 1, and 2, respectively. The term 5£eqcd in 
(0) includes all other local gauge-invariant operators of dimension 3 and higher that can 
be constructed out of A and A{. Static gauge-invariant correlation functions in full QCD 
can be reproduced in EQCD by tuning the gauge coupling constant g E , the mass parameter 
m E , the coupling constant X E , and the parameters in 5£eqcd as functions of g, T and the 
ultraviolet cutoff of EQCD. The Ag-dependence of the parameters is cancelled by the 
Ag-dependence of the loop integrals in the effective theory. 

In order to calculate the free energy using EQCD, we must also tune the coefficient 
of the unit operator, which was omitted from the effective lagrangian ([]) but appears as the 
first term in the expression (0) for the free energy. It depends on the ultraviolet cutoff A^ of 
EQCD in such a way as to cancel the cutoff dependence of the partition function for EQCD. 
The coefficient gives the contribution to the free energy from the momentum scale T. 
The logarithm of the partition function for EQCD includes the remaining contributions from 
the smaller momentum scales gT and g 2 T. 

In order to further separate the contributions from the scales gT and g 2 T, it is conve- 
nient to construct a second effective field theory called magnetostatic QCD (MQCD) which 
contains only the magnetostatic gauge field Af (x). The free energy for thermal QCD can be 
written 

F = T (f E (A E ) + f M (A E ,A M ) - ^QCD^) ; (1Q) 
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where ^mqcd is the partition function for MQCD: 

'(Am) 



^mqcd = J £>A?(x)exp(- J <2 3 x£ M qcd)- (11) 
The functional integral requires an ultraviolet cutoff Am- The lagrangian for MQCD is 



£mqcd — "opIjGlj + <5£mqcd, (12) 



where G?- is the magnetostatic field strength with coupling constant gu- This coupling 
constant differs from g E by perturbative corrections. The term <5£mqcd includes all possible 
local gauge-invariant operators of dimension 5 and higher that can be constructed out of Af. 
Gauge-invariant correlation functions in EQCD can be reproduced in MQCD by tuning the 
gauge coupling constant gu and the parameters in ££mqcd as functions of the parameters 
of EQCD (gE, m E , Xe, ■ ■ •) and the ultraviolet cutoff Am of MQCD. The A^-dependence 
of the parameters in the MQCD lagrangian is cancelled by the AM-dependence of the loop 
integrals in MQCD. 

In order to calculate the free energy using MQCD, one must also tune the coefficient fu 
of the unit operator, which was omitted from the effective lagrangian (|T2|) but appears as the 
second term in the expression fllPl) for the free energy. Its dependence on the ultraviolet cutoff 
Am of MQCD is cancelled by the cutoff dependence of the partition function for MQCD. 
The coefficient Jm gives the contribution to the free energy from the momentum scale gT. 
The contribution from the smaller momentum scale g 2 T is contained in the logarithm of the 
partition function for MQCD. 

By constructing the effective field theories EQCD and MQCD, we have separated the 
contributions from the momentum scales T, gT, and g 2 T in the free energy. The general 
structure of the free energy is 



T 



f E (T,g;A E ) + /j<f(m|, g E , As, . . . ; A E , A M ) + fa{9M, ■ ■ ■ ; Am) T, (13) 



where fa = —logZMQCT>/V . The arbitrary factorization scales A^ and Am separate the 
momentum scales T from gT and gT from g 2 T, respectively. The term f E and the parameters 



of EQCD (i.e., m 2 El g E , X E , . . .) involve only the scale T. They can therefore be calculated 
using ordinary perturbation theory as power series in g 2 (27cT), where g(2irT) is the running 
coupling constant at the scale of the lowest Matsubara frequency 2ttT. The term ju and 
the parameters of MQCD {gu-, ■ • •) involve only the scale gT. They can be calculated in 
EQCD as perturbation expansions in g 2 E /m E) X E /m E , and other dimensionless parameters 
obtained by multiplying EQCD coupling constants by appropriate powers of m E . The leading 
contribution to ju is proportional to m E . The term fa in ( p"3D can only be calculated using 
nonperturbative methods, such as lattice gauge theory simulations of MQCD. Surprisingly, 
however, Jq can be expanded as a weak coupling expansion in powers of g by treating the 
higher dimension operators in the MQCD lagrangian as perturbations J7|. The leading term 
is proportional to g\j. 

In summary, the free energy for QCD has the general structure given in (P^|). The term 
f E is the contribution from the scale T. It has the form T 3 multiplied by a power series in 
g 2 (27iT) whose coefficients can be calculated using ordinary perturbation theory in thermal 
QCD. The term J'm is the contribution from the scale gT. It has the form m E multiplied by 
a power series in g(2irT) whose coefficients can be calculated using perturbation theory in 
EQCD. The term fa is the contribution from the scale g 2 T. It has the form g\. { multiplied 
by a power series in g(27rT) whose coefficients can be calculated using lattice simulations of 
MQCD. 



3 Parameters in the EQCD Lagrangian 

In order to calculate the free energy using the EQCD lagrangian, the parameters in the 
lagrangian (^) must be tuned as functions of g, T, and so that EQCD reproduces the 
static gauge-invariant correlation functions of full QCD at distances R ^> 1/T. The EQCD 
parameters can be determined by computing various static quantities in full QCD, comput- 
ing the corresponding quantities in EQCD, and demanding that they match. It is convenient 
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to carry out these matching calculations using a strict perturbation expansion in g 2 . This 
expansion is afflicted with infrared divergences. The divergences arise from long-range forces 
mediated by static gluons, which remain massless in the strict perturbation expansion. Phys- 
ically, these divergences are screened by plasma effects either at the scale gT in the case of 
electrostatic gluons or at the scale g 2 T in the case of magnetostatic gluons. The screening 
of electrostatic gluons can be taken into account by summing up infinite sets of higher-order 
diagrams in the perturbation expansion, but the screening of magnetostatic gluons can only 
be taken into account using nonperturbative methods. Fortunately, it is not necessary to 
treat the effects of screening in a physically correct way in order to determine the parameters 
in the EQCD lagrangian. The parameters take into account the effects of large momenta of 
order T, and they are therefore insensitive to the infrared effects associated with screening. 
We can therefore simply ignore screening and remove the infrared divergences in the strict 
perturbation expansion by imposing any convenient infrared cutoff. As long as we use the 
same infrared cutoff in EQCD and in full QCD, we can determine the EQCD parameters 
by matching strict perturbation expansions in the two theories. Note that we are using 
the strict perturbation expansion simply as a device for determining the parameters in the 
EQCD lagrangian. 

3.1 Gauge coupling constant 

For the calculation of the free energy to order g 5 , we require the EQCD gauge coupling 
constant only to leading order in g 2 . At this order, we can simply read off ^from the 
lagrangian of the full theory. We substitute A (x,r) — > \^TA (x) in the QCD lagrangian 
(0) and compare Jq dr £qcd with £eqcd m (§)■ We find that, to leading order in g 2 , 

9 2 e = 9 2 T. (14) 

There is no dependence on the factorization scale at this order. The coupling constant 
qe could be calculated to higher order in g 2 by matching scattering amplitudes in full QCD 
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with the corresponding ones in EQCD. 



3.2 Mass parameter 

In this subsection, we calculate the coefficient m\ of the ^4q^o term in the EQCD lagrangian 
to next-to-leading order in g 2 . The physical interpretation of tue is that it is the contribution 
to the electric screening mass m e \ from large momenta of order T. The parameter m E can be 
determined by matching the strict perturbation expansions for the electric screening mass 
in full QCD and in EQCD. Beyond leading order in g, the electric screening mass becomes 



sensitive to magnetostatic screening and requires a nonperturbative definition [nj • However, 
in the presence of an infrared cutoff, m e i can be defined in full QCD by the condition that 
the propagator for the field Aq(t, x) at spacelike momentum K = (ko = 0, k) has a pole at 
k 2 = — m\. It is the solution to the equation 

k 2 + U{k 2 ) = at k 2 = -m 2 , (15) 

where U(k 2 ) is the \l — v — component of the gluon self-energy tensor evaluated at k = 0: 
Tl^ko = 0, k) = Tl(k 2 )5 ab . In EQCD with an infrared cutoff, the electric screening mass m e \ 
gives the location of the pole in the propagator for the field Aq(x). Denoting the self-energy 
function by UE(k 2 )5 ab , m el is the solution to 

k 2 + m 2 E + Ii E {k 2 ) = at k 2 = -m\. (16) 

By matching the expressions for m e i obtained by solving ( |T5] ) and (P!B[), we can determine 
the parameter m\. 

We calculate the electric mass m e \ in the full theory using a strict perturbation expansion 
in g 2 and using dimensional regularization with 3 — 2e spacial dimensions to cut off both 
infrared and ultraviolet divergences. The self-energy function U(k 2 ) can be expanded in a 
loop expansion 

U(k 2 ) = n (1) (A; 2 ) + n (2) (A; 2 ) + ... , (17) 
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with U^(k 2 ) and U^(k 2 ) being given by the diagrams in Fig. 1 and 2, respectively. We can 
simplify the equation (|15|) by expanding U(k 2 ) as a Taylor expansion around k 2 = 0. This is 
justified by the fact that the leading order solution to (|T5| ) gives a value of k 2 that is of order 
g 2 T 2 . The deviation of k 2 ^from should therefore be treated as a perturbation in order to 
get the strict perturbation expansion for m 2 ^ in powers of g 2 . The resulting expression for 
the electric screening mass to next-to-leading order in g 2 is 



m 



cl 



n«(o) + n( 2 )(o) - n«(o)^(o) . 



Here and below, we use the symbol "~" to denote an equality that holds only in the strict 
perturbation expansion. The one-loop diagrams that contribute to n^^A; 2 ) are shown in 
Fig. 1. Evaluating this function and its first derivative at k 2 = in Feynman gauge, we 
obtain 



II«(0) « Z g V|2(l-e)C A (J 1 -2J 1 ) - 4T F (J 1 -2J 1 )| 

2(l-e)(l + 2e 



(19) 



9 Z <-2C A 
+ 2T F 



X 2 + 



3-2e 



-J2 



3(1 -e) 
3-2e 



/C, 



~ 4(1 + 26) ~ 16 ~ 

J-2 H ^ 7: Jl — 7; TT^ 



3-2e 



3-2e 



(20) 



The sum-integrals X n , J7„, /C n X n , J/^, and /C n are defined in appendix A. The renormalization 
of the coupling constant using the MS scheme is accomplished by substituting 

l\C A -AT F g 2 1 



1 - 



(21) 



3 (An) 2 e 

into the expression for 11^(0). The two-loop diagrams that contribute to U(k 2 ) are shown 
in Fig. 2. This function evaluated at k 2 = is 

lK 2 )(0) « ( 7 4 |4(l- e )Cl[-2J 2l 7 1 + e J 1 X 2 + 4(l- e )J ll 7 2 ] 
+ 8C A T F [2X 2 J x - tZ{I 2 - 4(1 - e%J 2 

+ 8(1 - e)C F T F (Ti - Jx) (J 2 - 4J 2 ) |. (22) 
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The sum-integrals in (|T9|), (pO|), and fl22|) can be evaluated analytically using methods de- 
veloped by Arnold and Zhai [[HJ, and they are given in Appendix A. The three quantities 
appearing in (|T8| ) reduce to 



n«(o) 



n( 2 )(o) 



+ -7> 
3 



(An) 2 { 



1 + 

5 



"l + f 2 ^4 + 21og^ 
V C(-l) g 4vrT 

1-2^2 + 2 ,) ^ +21og 



c(-i; 



4vrT 



(23) 



C 



.4 



1 1 A 
--- + 2 7 + 21og-— 
e 5 Atti 



i-l+41og2 + 2 7 + 21og— - 
e Airl 



(24) 
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(4tt) 



T 2 



-C 2 



j +1+27+2 lriy +4l0g 4^ 



1 C'(-l) , A 

7 + 2-21og2 + 2 7 + 2^_i + 41og^ 



- ICpTp ,(25) 



where 7 is Euler's constant, C(z) is the Riemann zeta function, and A is the scale of dimen- 
sional regularization. Inserting these expressions into (|TBJ), we find that the strict perturba- 
tion expansion for m\ to order g A is 



rn 



el 



-g 2 (A)T 2 \C A + T F 



+ e 
+ 



22 



4ttT, 
22 , A 



+ Tp l-21og2 + 2 



16 



^U + y^y 10 ^; + ^ T H 3 -y log2+ y 7+ y log 4^Tj 



C'(-i) 
C(-i) 

14 



+ 2 log 



14 



A 
4ttT 



A \ 



A \ 



4 16 8 8 

log 2 7 log 

3 3 & 3 ' 3 & 4vrT 



J 



6C F T F 



5 



(47T) 



(26) 



Note that all the poles in e have cancelled. In the order g 2 term, we have kept terms of order 
e for later use. The expression (|26"D depends on A explicitly through logarithms of A/AirT 
and implicitly through the coupling constant <? 2 (A). The scale of the coupling constant 
can be shifted from the dimensional regularization scale A to an arbitrary renormalization 
scale /1 by using the solution to the renormalization group equation for the running coupling 
constant: 



2 (A) = g 2 M 



2(11C A -4T F ) g 2 ii 
1 + 3 (4^ l0g A 



(27) 
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After making this shift in the scale of the coupling constant, the only remaining dependence 
on A occurs in the terms of order e. In these terms, A can be identified with the factorization 
scale that separates the scales T and gT. 

The expression (EJ]) for m\ is an expansion in powers of g 2 . It does not include a g 3 
term, in contrast to the expression for m\ that correctly incorporates the effects of the 
screening of electrostatic gluons [ES]. This g 3 term arises because the g 4 correction includes 



a linear infrared divergence that is cut off at the scale gT. Since we have used dimensional 
regularization as an infrared cutoff, power infrared divergences such as this linear divergence 
have been set equal to 0. 

In order to match with the expression (p6|), we have to calculate the screening mass in 
EQCD using the strict expansion in g 2 . Since m\ is treated as a perturbation parameter of 
order g 2 , the only scale in the self-energy function n^/c 2 ) is k 2 . After Taylor expanding in 
powers of k 2 , there is no scale in the dimensionally regularized integrals, so they all vanish. 
The solution to the equation flTED for the screening mass is therefore trivial: 

m\ m|. (28) 



Comparing (|26|) and (28) , we find that, in the limit e — > 0, the parameter m 2 E is given by 
1 



m 



E 



e=0 

+ 



3 S 2 (/i) T 2 lC A + T F 



^ 2 /5 22 22 . y. \ _ _ / 16 . n 14 14 , p \ 

c * (3 + T 7 + T hg £f) + CaTf ( 3 " T log2 + T 7 + T log i^) 



+ ^l--^log2-^-^log/-l - QCpTp 



I 1(5 , , s s , ^ \ ^] g 2 



(4tt) 



(29) 



.3 3 3' 3 °47rTy 
At this order in g 2 , there is no dependence on the factorization scale A#. The order-e terms 
in m\ will also be required later in the calculation. These terms are 



dm% 



de 



3" [ v C(-i) 4ttT; 

+ T F (l-21og2 + 2^ + 21og^U (30) 



This expression depends explicitly on the factorization scale A E . 

14 



3.3 Coefficient of the unit operator 



In this subsection, we calculate the coefficient of the unit operator to next-to-next-to- 
leading order in g 2 . The physical interpretation of is that /#T is the contribution to the 
free energy from large momenta of order T. The parameter f F is determined by calculating 
the free energy as a strict perturbation in g 2 in both full QCD and EQCD, and matching 
the two results. 

In the full theory, the free energy has a diagrammatic expansion that begins with the one- 
loop, two-loop and three-loop diagrams shown in Fig. 3, 4, and 5. Evaluating the diagrams 
in Feynman gauge, we obtain 

F « -(1 - e)d A T + 2d F l' 

+ d A Z 2 g 2 \{l - e) 2 C A l\ + 2(1 - e)T F l x {l x - 21,) 



+ d A g^C 2 A (l - e) 2 
+ C A T F (1 - e) 



2(1 + e)l 2 l 2 - -M ,o ~ M 2 ,- 2 



-8JiX 2 Xi -2eM 0t 



e)Wo, + AM -2,2 



8(1 + e)l 2 l 2 + 2eA/"o, - 4A/" 2 ,_ 2 



2C F T F (l - e) 2(1 - e) (if - 4T 1 T 1 + if) 1 2 + 2M , - (1 + e)A/" , 



+2(l-e)M lr 



(31) 



The symbol "~" is a reminder that the strict perturbation in g 2 does not give a physically 
correct treatment of the screening effects of the plasma. The sum-integrals in (|31~1) are given 
in Appendix |A|. To order g 4 , the renormalization of the coupling constant is accomplished 
in the MS scheme by substituting (|2lD for Z g in the order g 2 term. The final result is 



ii 2 d ArTl4 \l , 7 d f („ , 5^ ^ g 2 {A) 



+ 
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In EQCD, the free energy is given by the expression (|7|). We calculate log Zeqcd using the 
strict perturbation expansion in which g 2 E and m 2 E are treated as perturbation parameters and 
both infrared and ultraviolet divergences are regularized using dimensional regularization. 
Since diagrams with massless propagators and with no external legs vanish in dimensional 
regularization, the only contribution to logi?EQCD which does not vanish comes from the 
counterterm 5f E which cancels ultraviolet divergences proportional to the unit operator. 
The resulting expression for the free energy is simply 

F ~ (/# + Sf E ) T . (33) 

The counterterm can be determined by calculating the ultraviolet divergences in logZ E qcD- 
If we use dimensional regularization together with a minimal subtraction renormalization 
scheme in the effective theory, then SJe is a polynomial in g E , m E , and the other parameters 
in the lagrangian for EQCD. The only combination of parameters that has dimension 3 and 
is of order g 4 is g E m E . Thus the leading term in Sf E is proportional to g E m 2 E . The coefficient 
is determined by a 2-loop calculation that is a trivial part of the 3-loop calculation in Section 
4. The result for the counterterm is 

6fE = ~M4^ gEmE ~e • (34) 

When expressing this counterterm in terms of the parameters g and T of the full theory, 
we must take into account the fact that m 2 E multiplies a pole in e. Thus in addition to 
expression for m 2 E given in fl29|) , we must also include the terms of order e which are given 
by (j30|). The counterterm (|34]) is therefore 

12Ci(l + 2M + 21ogM 



5h - ~ m~ 2 



C(-l) °4ttT / 



+ \2C A T F ( - + 1 - 2 log2 + 2^-^ + 2 log Ai 
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Note that minimal subtraction in the effective theory is not equivalent to minimal subtraction 
in the full theory. In addition to the poles in e in (|35|), there are finite terms that depend on 
the factorization scale A F . 

Matching (|32|) with ( p3]) and using the expression (S3J), we conclude that to order g A 



is 



f E (A E ) 



7t 2 d 
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48 log 
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(36) 



where g(/i) is the coupling constant in the MS renormalization scheme at the scale /i. We 
have used ( f27j ) to shift the scale of the running coupling constant from A^ to an arbitrary 
renormalization scale fi, and we have identified the explicit factors of A that remain with 
the factorization scale A F - 



3.4 Evolution of EQCD coupling constants 

The effective lagrangian (Q) for EQCD can be expressed as a sum over all local operators 
that respect the symmetries of the theory: 



f E (A E ) + £ EQ cd = J2 C n(A E )O n , 



(37) 



where we have included the unit operator with coefficient f E as one of the operators O n . 
The coefficients C n are the generalized coupling constants of the effective theory. Because 
of ultraviolet divergences, the effective theory must be regularized with an ultraviolet cutoff 
A^. The ultraviolet divergences in the effective theory include power ultraviolet divergences 
proportional to A^, p = 1, 2, . . ., and logarithmic divergences proportional to log(A£). The 
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power divergences are artifacts of the regularization scheme and have no physical content. 
If they are not removed as part of the regularization procedure, they must be cancelled 
by power divergences in the coupling constants C n . In contrast, the logarithmic ultraviolet 
divergences are directly related to logarithms of T in the full theory, and therefore represent 
real physical effects. It is convenient to use a regularization procedure for the effective theory 
in which power ultraviolet divergences are automatically subtracted, such as dimensional 
regularization. In this case, they need not be cancelled by power divergences in the coupling 
constants. The dimensions of a coupling constant can then only be taken up by powers of the 
temperature T. The coupling constant C n must be proportional to T 3 ~ dn , where d n is the 
scaling dimension of the corresponding operator O n . The dimensionless factor multiplying 
rpz-dn j n COU pii n g constant C n can be computed as a perturbation series in g 2 {T), with 
coefficients that are polynomials in \og(T / h E ). The dependence on is governed by a 
"renormalization group equation" or "evolution equation" of the form 

A E ^-C n (A E ) = (3 n (C), (38) 

where the beta function j3 n has a power series expansion in the coupling constants C m . 
These equations follow from the condition that physical quantities must be independent of 
the arbitrary scale A E . Since C n is proportional to T 3 ~ dn , every term in the expansion of its 
beta function must be proportional to T 3 ~ dn . In particular, a term like C mi C m2 . . . C mk can 
appear only if the dimensions d m . of the corresponding operators O mi satisfy 

k 

£(3-d m J = 3-d n . (39) 

8=1 

The condition (|39"D is very restrictive, particularly if the effective lagrangian is truncated 
to the super- renormalizable terms that are given explicitly in (|J). It implies that the only 
terms that can appear in the beta function for the coefficient f E of the unit operator are 
g E m E , X E m E , and a cubic polynomial in g E and X E . Since m E , g%, and X E are of order g 2 , g 2 , 
and g 4 , respectively, the only term of order g A is g 2 E m 2 E . We can determine its coefficient by 
calculating the ultraviolet divergences in the strict perturbation expansion for the free energy 
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in the effective theory. These divergences do not appear in (p3|), because the ultraviolet poles 
in e have cancelled against infrared poles in e. We can calculate the ultraviolet divergences 
by using a different regularization for infrared divergences. Alternatively, since we have 
already calculated f E explicitly to order g 4 , we can simply differentiate (|36| ) and use the fact 
that A E (d/dA E )f E must be proportional to g 2 E m 2 E . Using g 2 E = g 2 T and the leading order 
expression for m% in (p9"D, we find that the evolution equation is 

Ae^-Je = ~^9We + O(g^) . (40) 

The beta function for m E must be a quadratic polynomial in g\ and X E . The terms 
g E , g E X E , and X 2 E are of order g A , g 6 , and g 8 , respectively. The coefficients of these terms 
can be determined by calculating the ultraviolet divergent terms in the strict perturbation 
expansion for the electric screening mass in the effective theory. Alternatively, if m 2 E is 
known, its beta function can be determined simply by differentiating. Since the expression 
([29]) is independent of A E , we know that the coefficient of g% in the beta function vanishes 
and the leading term must be g 2 E \ E . Thus the evolution equation for m 2 E is 

A^m| = O + O^T 2 ). (41) 

We have not calculated the coefficient of g E \ E in this evolution equation, because it does 
not affect the free energy until order g 7 . 

The beta functions for g 2 E and vanish to all orders in the super-renormalizable inter- 
actions. All the nonvanishing terms in their beta functions involve the coupling constants of 
nonrenormalizable interactions, and they are therefore suppressed by large powers of g. The 
evolution of these paramaters can probably be ignored for most practical purposes. 

The only EQCD parameter whose evolution plays a role in the free energy to order g 6 is 
f E . To this order, the solution to the equation (|40| ) is trivial: 



f E (A E ) = f E (A' E ) - ^gWE log ^ • (42) 
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4 Free Energy to Order g 5 

Having calculated the parameters of EQCD to the necessary order in g 2 , we now use the 
effective theory to calculate the free energy to order g 5 . The free energy is the sum of the 
three terms in (Jl3|), which correspond to the momentum scales T, gT, and g 2 T, respectively. 
The term JeT is the contribution from the scale T. We have already calculated Je to the 
necessary order and it is given in fl36|). The term faT is the contribution from the scale g 2 T, 
but it does not contribute until order g 6 . The remaining term JmT is the contribution ^from 
the scale gT. 

Through order g 5 , ju is proportional to the logarithm of the partition function for EQCD: 

log ^EQCD /.oN 

JM y • (4dJ 

In order to calculate this contribution using perturbation theory, we must incorporate the 
terms in the lagrangian that provide electrostatic screening into the free part of the la- 
grangian. The necessary screening effects are provided by the AqAq term in the EQCD 
lagrangian. Thus we must include the effects of the mass parameter m 2 E to all orders, while 
treating all the other coupling constants of EQCD as perturbation parameters. The only 
coupling constant that is required to obtain the free energy to order g 5 is the gauge coupling 
constant g^. 

The contributions to logZ E Q CD of orders g 3 , g 4 , and g 5 are given by the sum of the 1- 
loop, 2-loop, and 3-loop diagrams in Fig. 6, 7, and 8. The solid, wavy, dashed lines represent 
the propagators of the A field, the Aj fields, and the associated ghosts, respectively. We 
evaluate these diagrams in Feynman gauge. They can be expressed in terms of the scalar 
integrals defined in Appendix B. The resulting expression for the logarithm of the partition 
function is 



fu = ~^~Iq + ^aCa9e 



\ l l + m 2 E Ji 
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(44) 



where is the counterterm associated with the unit operator of the EQCD lagrangian. 
The integrals I n , J n , K n , L m>n , M m ^ n , and N m ^ n can be calculated analytically using methods 
developed by Broadhurst [23] and they are given in Appendix B. Adding them up, we obtain 
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where A is the scale of dimensional regularization. It can be identified with the ultraviolet 
cutoff of EQCD. The ultraviolet pole in e in the term proportional to g%m E in J45j) is 



cancelled by the counterterm Sf E , which is given in (34j). Our final result is therefore 
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The coefficient /m in (f£E| ) can be expanded in powers of g by setting gf|, = g 2 T and 
by substituting the expression ( p9|) for m|;. The complete free energy to order g 5 is then 
F — (/b + fhi)T. Note that the dependence on the arbitrary factorization scale A s cancels 
between and Jm, up to corrections that are higher order in g, leaving a logarithm oiT/m E . 
This g 4 log(g) term is associated with the renormalization of f E , and its coefficient can be 
determined from the evolution equation (|4"0|). There is no g 5 log(g) term in the perturbation 
expansion for F, and this is a consequence of the vanishing of the order-g 4 term in the beta 



function for m% 
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5 Outline of Calculation to Order g 6 

The calculation of the free energy to order g 5 , which was presented in the previous section, 
was greatly streamlined by using effective field theory to unravel the effects of the momentum 
scales T and gT. The same result has also been obtained by Kastening and Zhai using other 
methods ]I8| . However the advantages of the effective field theory approach become more 
and more apparent as we go to higher order in g. In this section, we demonstrate the power 
of this method by outlining the calculation of the free energy to order g G . In this case there 
are contributions from all three momentum scales T, gT, and g 2 T. 

5.1 Contribution from the scale g 2 T 

We first discuss the contribution to the free energy from the scale g 2 T, which is given by the 
term f@T in (0). This term is proportional to the logarithm of the partition function ( |i~l~D 
of MQCD. Treating the correction term 5£mqcd in the MQCD lagrangian as a perturbation, 
the partition function can be written 

^mqcd = f AM ' W(x)e X p(- J d 3 xG 2 /4){l - Jd 3 x5C MQC B + •••}, (47) 
where G 2 = G^G^. Taking the logarithm of both sides, we obtain 

fo = -^f^ + (5C MQCD ) + . . . , (48) 

where ^mqcd * s the partition function for the minimal gauge theory with action / d?xG 2 /A. 
The subscript on the expectation value (<5£mqcd)o is a reminder that it is to be calculated 
using the minimal gauge theory action. 

For the moment, let us consider only the first term in fl48l) . The partition function 
^mqcd is that of the minimal gauge theory in 3 dimensions. This is a super-renormalizable 
theory and its ultraviolet divergences have a very simple structure. By naive power-counting, 
ultraviolet divergences in log^^Q CD can arise only from vacuum diagrams with 1, 2, 3, or 
4 loops or from propagator corrections with 1 or 2 loops. Ward identities guarantee that 

22 



the propagator corrections are actually finite. This is related to the fact that the only 
gauge invariant operator with dimension lower than G 2 is the unit operator. Thus the 
only ultraviolet divergences are in the vacuum diagrams. The 1-loop diagrams give a cubic 
divergence. The 2-loop diagrams give a quadratic divergence proportional to g M . The 3-loop 
diagrams give a linear divergence proportional to g A M . Finally, the 4-loop diagrams give a 
logarithmic divergence proportional to g M . After subtraction of the power divergences, we 
can use dimensional analysis to determine the form of log Zmqqd • Aside from the logarithmic 
dependence on the ultraviolet cutoff Am, the only scale in the problem is gM- By dimensional 
analysis, logZ^Q CD must be proportional to g M . Thus it must have the form 



where a and b are pure numbers. The coefficient b can be determined by calculating the 
logarithmic ultraviolet divergence in the 4-loop vacuum diagrams for MQCD. The coefficient 
a can only be calculated using nonperturbative methods. It can for example be extracted 
from measurements of the expectation value (G 2 )o using lattice simulations of the pure gauge 
theory. A convenient expression for (G 2 )o can be obtained by taking the logarithm of both 
sides of ( |4lD and differentiating with respect to g 2 M . It is useful to first rescale the field 
Ai in the functional integral for ^mq CD , so that the coupling constant appears only in the 
coefficient l/g M of the action. After subtracting the power ultraviolet divergences, we obtain 
the expression 



The subscript on the expectation value (G 2 )q is a reminder that it is to be calculated 
using the minimal gauge theory action J d 3 xG 2 /4 rather that the full action of MQCD. The 
expectation value (G 2 )q can be measured on the lattice using Monte Carlo simulations of 
the minimal gauge theory. Once (G 2 )q has been measured and the coefficient b has been 
calculated, we can determine a using the formula (|50|) . 

We now verify that the correction term in ([17]) from higher dimension operators in the 




(49) 




(50) 
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MQCD lagrangian can indeed be treated as a small perturbation. The lowest dimension 
operators in <5£mqcd are G 3 = f ahc G\jG h ^ k G c ki ^ whose coefficient is proportional to g 3 /T 3 ^ 2 , 
and (DG) 2 = (DiGik) a {DjGjk) a , whose coefficient is proportional to g 2 /T 2 . Their coefficients 
have been calculated to leading order in g by Chapman for the case of a pure gauge theory 
pOfl . After subtraction of power ultraviolet divergences, the only scale in the problem is g 2 M . 
Therefore, by dimensional analysis, (G 3 ) must be proportional to (<?m) 9//2 . Using g 2 M m g 2 T 
and taking into account the coefficient which is proportional to g 3 /T 3 ^ 2 , we find that the 
contribution to fa from (G 3 )o is of order g 12 T 3 . Using a similar analysis, we find that 
the contribution from {(DG) 2 )q is also of order g 12 T 3 . Thus the effects of higher dimension 
operators in the MQCD lagrangian are indeed suppressed by powers of the coupling constant 
9- 

We have found that the contribution to the free energy from the scale g 2 T can be written 

foT = (a + 6 log ^] g\[T + 0(s 12 T 4 ) , (51) 
V 9 m J 

Remarkably, the only nonperturbative calculation that is required to determine the free 
energy up to order g 12 is that of the single pure number a. We also require the coupling 
constant gM, which can be calculated by matching perturbative calculations in EQCD and 
MQCD. To calculate the free energy to order g 6 , we only need gM to leading order in g. At 
this order, it is given simply by g 2 M = g 2 T. In summary, in order to obtain the contribution 
to the free energy from the scale g 2 T to order g 6 , all that is required are the two pure numbers 
a and b in (|5T|). The number b can be calculated by evaluating 4- loop diagrams in MQCD. 
In Ref. J7J, it was assumed incorrectly that this number vanishes. The number a can be 
calculated using lattice simulations of the pure gauge theory in 3 dimensions. 

5.2 Contribution from the scale gT 

The contribution to the free energy from the scale gT is given by the term JmT in (fl~3|). 
The coefficient Jm can be determined by calculating the logarithm of the EQCD partition 
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function in both EQCD and MQCD and matching the expressions. If we use dimensional 
regularization to cut off both infrared and ultraviolet divergences, all the loop diagrams in 
MQCD vanish. The expression for the logarithm of the partition function then is simply 

lQgZ E QCD , Af ft .„. 
y - JM + OjM , (,52) 

where 5/m is a counterterm that cancels ultraviolet divergences in MQCD that are propor- 
tional to the unit operator. To order g 6 , this counterterm is simply 

SfM = y/m , (53) 

where b is the same coefficient that appears in (f)Tl). 

To determine fu-, we must match the expression ( |52|) with the corresponding expression 
in EQCD, which is obtained by calculating the sum of vacuum diagrams using dimensional 
regularization. The resulting expression for logZ E Q CD is a sum of polynomials in the EQCD 
coupling constants, such as g\ and Xe, multiplied by whatever powers of mg are required 
by dimensional analysis. There are three such terms that contribute to the free energy at 
order g 6 . The first term is g E m E , whose coefficient has already been calculated in fl46|) . 
It contributes through the next-to-leading order term in m E , which is given in (|29|) , and 
through the next-to- leading order term in g E , which has not yet been calculated. The 
second term which contributes at order g 6 is proportional to g E . Its coefficient is determined 
by calculating all 4- loop vacuum diagrams that involve only the gauge coupling constant gE- 
This term will have a pole in e that matches that from the counterterm ([53]). The third 
term that contributes to ju at order g 6 is proportional to A^m^. Its coefficient is given by 
the single 2-loop vacuum diagram that involves the Aq coupling constant Xe- This coupling 
constant is only required to leading order in g and has already been calculated by Nadkarni 
and by Landsman |2l], |J. 

In summary, there are three coefficients that must be calculated in order to obtain the 
contribution of order g 6 to the free energy from the scale gT. We need the coefficients of g% 
and of Xettie in the expression for fu- These can be obtained by perturbative calculations 
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in EQCD. We also need the coefficient of g A in the expression for the EQCD parameter g\. 
This requires a perturbative calculation in full QCD. 



5.3 Contribution from the Scale T 



The contribution to the free energy from the scale T is given by the term /#T in (13). The 
term is obtained by matching the strict perturbation expansions for the free energy in 
full QCD and in EQCD. In full QCD, the contribution of order g 6 is the sum of all 4-loop 
vacuum diagrams. If we use dimensional regularization to cut off both infrared and ultraviolet 
divergences, then the corresponding expression in EQCD is simply F = (/# + 5/b)T. The 
counterterm 8/e includes the term proportional to g 2 E m 2 E /e given in ( [34]) and also a term 
proportional to Agml/e. Since the counterterm is proportional to 1/e, we need not only the 
value of the coupling constant Xe at e = but also the terms linear in e. Similarly, we need 
the term of order e in the order- g 4 correction to g\. 

In summary, there are several calculations that must be carried out in order to obtain 
the term of order g 6 in f&. We need to calculate the 4-loop vacuum diagrams in full QCD. 
We also need to calculate the terms of order eg 4 in the EQCD parameters g\ and Xe- 



6 Convergence of Perturbation Theory 

We have calculated the free energy as a perturbation expansion in powers of g to order g 5 . 
In this section, we examine the convergence of that perturbation expansion. For simplicity, 
we focus on the case of QCD with n/ flavors of quarks. 

The effects of the momentum scale T enter into the free energy only through the coefficient 
Je and the parameters in the EQCD lagrangian. The term is given in (j36p: 
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+ 0{al) } . (54) 



The other parameters in the EQCD lagrangian that enter into the calculation of the free 
energy to order g 5 are rriE and gE, which are given by (p9|) and (|T4]), respectively: 



m, 



4vr a s (ji) T 2 |l + \n f + 0.612 - 0.488^ - 0.0428 
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+ 0(« s 2 ) , (55) 



g% = 4vr a s T [1 + 0(« s )] . (56) 

We have calculated two terms in the perturbation series for m 2 E and three terms in the 
series for f E . We can use these results to study the convergence of perturbation theory for 
the parameters of EQCD. We consider the case of n/ = 3 flavors of quarks, although our 
conclusions will not depend sensitively on rif. The question of the convergence is complicated 
by the presence of the arbitrary renormalization and factorization scales \x and Ae- The next- 
to- leading-order (NLO) correction to Je is independent of \x and Ae, and is small compared to 
the leading-order (LO) term provided that a s (fi) < 1.1. The NLO correction to m E and the 
next-to-next-to-leading-order (NNLO) correction to J'e both depend on the renormalization 
scale jji,. One scale-setting scheme that is physically well-motivated is the BLM prescription 
||26|| , in which fi is adjusted to cancel the highest power of rif in the correction term. This 
prescription gives /z = 0.93 ttT when applied to m 2 E and fi = 4AnT when applied to Je- 
These values differ only by about a factor of 2 from 2ttT, which is the lowest Matsubara 
frequency for gluons. Below, we will consider the three values fi = ttT, 2ttT, and 4ttT '. For 
the NLO correction to m 2 E to be much smaller than the LO term, we must have a s (fi) 
0.8, 3.8, and 1.4 if \i = nT, 2tiT, and 4nT, respectively. Based on these results, we conclude 
that the perturbation series for the parameters of EQCD are well-behaved provided that 
a s (2uT) < 1. 

The NNLO correction for f E depends not only on /i, but also on the factorization scale 
A^. Because the coefficient of log(Ae/27rT) in (|36|) is so much larger than that of log(/z/27rT), 
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the NNLO correction for f E is much more sensitive to than to /i. It is useful intuitively to 
think of the infrared cutoff Ag as being much smaller than the ultraviolet cutoff /i. However, 
these scales can be identified with momentum cutoffs only up to multiplicative constants 
that may be different for /x and A#. Both parameters are introduced through dimensional 
regularization, but \x arises from ultraviolet divergences of 4-dimensional integrals, while Ae 
arises ^from infrared divergences of 3-dimensional integrals. We might be tempted to set 

= [i, but then the NNLO coefficient in f E is large. For the choice \i = 2ttT, the correction 
to the LO term is a multiplicative factor 1 — 0.9a s + 6.5a;;. The NNLO correction can be 
made small by adjusting Ae- It vanishes for Ae = 5.8 irT, 5.1 ttT, and 4.5 7rT if \x = tiT, 
2nT, and AtiT, respectively. We conclude that the perturbation series for is well-behaved 
if the factorization scale A^ is chosen to be approximately 5ttT. Whether this choice is 
reasonable can only be determined by calculating other EQCD parameters to higher order 
to see if the same choice leads to well-behaved perturbation series. 

The choice of A^ that makes the perturbation series for the EQCD parameters well- 
behaved may be much larger than the largest mass scale uie of EQCD. Perturbative correc- 
tions in EQCD will then include large logarithms of Ae/tue- This problem can be avoided 
by using renormalization group equations to evolve the parameters of EQCD from the initial 
scale Ae down to some scale A' E of order m^. The solution to the renormalization group 
equation for is given in (f4*2"|). The evolution of g E and m 2 E occurs only at higher order in 
the coupling constant and therefore can be ignored. 

We have carried out only one perturbative calculation in EQCD. This is the term J'm, 
which gives the contribution to the free energy from the scale gT. This term is given in 

2 



2 

/m(A E ) = -7^ m E 



1 - (0.256-^log^)-^ - 27.6 + Otf 



(57) 

We now consider the convergence of the perturbation series fl46|) for Jm- The size of the NLO 
correction depends on the choice of the factorization scale Ae- It is small if A^ is chosen to 



be approximately mg. The NNLO correction in (f46D is independent of any arbitrary scales. 
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If n/ = 3, it is small compared to the leading order term only if a s <C 0.17. Thus the 
perturbation series for Jm is well-behaved only for values of a s (27rT) that are much smaller 
than those required for the parameters of EQCD to have well-behaved perturbation series. 

Inserting (|55|) and ([56]) into (|57|), expanding in powers of g, and adding (0), we get the 
expansion for the free energy in powers of JoT s : 
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The coefficient F 2 was first given by Shuryak . The coefficient of F 3 was first calculated 
correctly by Kapusta [23||. The coefficient F 4 was calculated in 1994 by Arnold and Zhai 



13|. The coefficient F 5 has also been calculated independently by Kastening and Zhai [18 



We now ask how small a s must be in order for the expansion (|58|) to be well-behaved. For 
simplicity, we consider the case nj = 3, although our conclusions are not sensitive to rif. If we 
choose the renormalization scale /i = 2ttT motivated by the BLM criterion | 26| , the correction 
to the LO result is a multiplicative factor 1 - 0.9a s + 3.3a^ /2 + (7.1 + 3.5 \oga s )a 2 s - 20.8c^ /2 . 
The oPj 2 term is the largest correction unless a s (27xT) < 0.12. We can make the ct'J 2 term 
small only by choosing the renormalization scale to be near the value [i = 36.5nT for which 
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F 5 vanishes. This ridiculously large of ji arises because the scale \i has been adjusted to 
cancel the large g 5 correction to J'm in (0)- This contribution arises from the momentum 
scale gT and has nothing to do with renormalization of a s . We conclude that the expansion 
(|58|) for F in powers of y^a^ is well-behaved only if a s (2irT) <C 1/10. This is an order of 
magnitude smaller than the value required for the EQCD parameters to be well-behaved. 
Our previous analysis indicates that this slow convergence of the expansion for F in powers 
of y/al can be attributed to the slow convergence of perturbation theory at the scale gT. 

7 Discussion 

In this paper, we have used effective-field-theory methods to unravel the contributions to 
the free energy of high temperature QCD from the scales T, gT, and g 2 T. We calculated 
the free energy explicitly to order g 5 . The calculation was significantly streamlined by using 
effective-field-theory methods to reduce every step of the calculation to one that involves 
only a single momentum scale. We also outlined the calculations that would be necessary 
to obtain the free energy to order g 6 . It is only at this order that the full power of the 
effective- field-theory approach becomes evident. 

The effective-field-theory approach provides an understanding of the logarithms of the 
coupling constant that arise in perturbation expansions in thermal field theory. These loga- 
rithms are associated with the renormalization of the parameters of effective field theories. 
The resulting evolution equations can be used to sum up leading logarithms of the coupling 
constant of the form g m+2n log n (g) to all orders in n ||. To the accuracy required for the 
calculation of the free energy to order g 5 in QCD, this resummation is trivial. The only 
terms of the form g m+2n \og n (g) with m + 2n < 5 are a g A log(g) term associated with renor- 
malization of the coefficient /g. The fact that the solution ( fE|) to the evolution equation 
for /e is trivial indicates that there are no higher order terms of the form g 2+2n log ra (g) that 
are related to the g A log(g) term through the renomalization group. There are also no terms 
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of the form g 3+2n log n (g) in the free energy. This is a consequence of the vanishing of the 
g E term in the beta function for m? E . In the seemingly simpler problem of a massless scalar 
field with a 4 interaction, the evolution equations play a more important role [Q. There are 
terms in the free energy of the form g 3+2n log n (g) that can be summed up to all orders with 
the help of the renormalization group. The relative simplicity of the QCD case comes from 
the fact that the term g E in the beta function for m 2 E has a vanishing coefficient. We know 
of no deep reason for this coefficient to vanish. 

Our explicit calculations allow us to study the convergence of the perturbation expansion 
for thermal QCD. They suggest that perturbation theory at the scale gT requires a much 
smaller value of the coupling constant than perturbation theory at the scale T. At the 
scale T, perturbation corrections can be small only if a s (2irT) < 1. Of course, even if 
this condition is satisfied, the perturbation expansion may break down anyway, but this 
is certainly a necessary condition. At the scale gT, perturbation corrections can be small 
only if a s (27rT) <C 1/10. Thus, in order to achieve a given relative accuracy, the coupling 
constant a s (2nT) must be an order of magnitude smaller for perturbation theory at the 
scale gT compared to perturbation theory at the scale T. This has important implications 
for calculations in thermal QCD. At extremely high temperatures, the asymptotic freedom 
of QCD guarantees that the running coupling constant a s {2nT) is sufficiently small that 
perturbation theory will provide an accurate treatment of the effects of the scale gT as well 
as those of the scale T. Nonperturbative methods, such as lattice simulations of MQCD, are 
necessary only to calculate the effects of the scale g 2 T. Of course, one can always treat the 
entire problem nonperturbatively by carrying out lattice simulations of full thermal QCD. 
However it is probably more efficient to integrate out the scales T and gT using perturbative 
methods, and to reserve the nonperturbative methods only for the scale g 2 T where they 
are essential. As the temperature is decreased, the running coupling constant increases and 
perturbation theory becomes less accurate. At sufficiently low temperatures, perturbation 
theory breaks down completely, and the entire problem must be treated nonperturbatively. 
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This is certainly the case when the temperature is close to the critical temperature for the 
phase transition from the quark-gluon plasma to a hadron gas. 

Our calculations suggest, however, that there is a range of temperatures in which per- 
turbation theory at the scale gT has broken down, but perturbation theory at the scale T 
is reasonably accurate. In this case, one can still use perturbation theory at the scale T to 
calculate the parameters in the EQCD lagrangian. Our calculations of the coefficients 
and m 2 E to order g A are therefore still useful. However, nonperturbative methods, such as 
lattice simulations of EQCD, are required to calculate the effects of the smaller momentum 
scales gT and g 2 T. While one could simply treat the entire problem nonperturbatively us- 
ing lattice simulations of full QCD, the effective-field-theory approach provides a dramatic 
savings in resources for numerical computation. The savings come ,Trom two sources. One 
is the reduction of the problem from a 4-dimensional field theory to a 3-dimensional field 
theory. The other source of savings is that quarks are integrated out of the theory, which 
reduces it to a purely bosonic problem. 

We now consider briefly the implications for the study of the quark-gluon plasma in 
heavy-ion collisions. The critical temperature T c for formation of a quark-gluon plasma is 
approximately 200 MeV. It may be possible in heavy-ion collisions to produce a quark-gluon 
plasma with temperatures several times T c . At T = 350 MeV, a s (27iT) pa 0.3, which is 
small enough that perturbation theory may be reasonably convergent at the scale T, but it 
is certainly not convergent at the scale gT. We conclude that at the temperatures achievable 
in heavy-ion collisions, perturbative QCD may be accurate when applied to quantities that 
involve the scale T only. However nonperturbative methods are required to accurately calcu- 
late quantities that involve the scales gT and g 2 T. The most effective strategy for calculating 
the properties of a quark-gluon plasma at such temperatures will probably involve a com- 
bination of perturbative and nonperturbative methods. The effective-field-theory approach 
developed in this paper provides a systematic method for unraveling the momentum scales 
in the plasma and for combining perturbative and nonperturbative methods in a consistent 
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way. This approach applies strictly only to static properties and to the case of zero baryon 
density. The extension to dynamical properties and to the case of nonzero baryon density 
remains a challenging problem. 
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A Sum- integrals in the Full Theory 

In the imaginary-time formalism for thermal field theory, the 4-momentum P = (po,p) is 
euclidean with P 2 = p^ + p 2 . The euclidean energy p has discrete values: po = 2niiT for 
bosons and po = (2n + l)nT for fermions, where n is an integer. Loop diagrams involve 
sums over po and integrals over p. It is convenient to use dimensional regularization to 
regularize both ultraviolet and infrared divergences. We introduce a concise notation for 
these regularized sum-integrals: 



where 3 — 2e is the dimension of space and A is an arbitrary momentum scale. The factor 
(e 7 /47r) e is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, A coincides with the renormalization scale in the MS renormalization scheme. 
Below, we collect together all the sum-integrals that are required to calculate the coefficient 
Je to next-to-next-to-leading order in g 2 and the coefficient m 2 E to next-to-leading order in 




(A.2) 



(A.l) 
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The one-loop bosonic sum-integrals that arise in the calculation have the following forms: 

1 



p (P 2 )™' 



pi 
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The specific sum-integrals that are needed are 



4 

Po 



(A.3) 
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2c 
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(A.6) 
(A.7) 

(A.8) 

(A.9) 

(A.10) 

(All) 



where 7 is Euler's constant and ((z) is Riemann's zeta function. In ( |A.6| ), I' denotes the 
derivative of X n with respect to n evaluated at n = 0. The one- loop fermionic sum-integrals 
have the following forms: 
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T 



{P} (P 2 )"' 



{P} (P 



2\n+l 

4 
PI 



t{P} (P 2 )"+ 2 

The specific sum-integrals that are needed are 

7tt 2 
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(A.13) 
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(A.16) 
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J ' = if 
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(A.17) 
(A.18) 
(A.19) 
(A.20) 



All of the two-loop sum-integrals that arise in the calculation factor into the product of 

2 one-loop sum-integrals. Some of the three-loop sum-integrals factor into the product of 

3 one-loop sum-integrals. Others factor into the product of a one-loop sum-integral and a 
two-loop sum-integral. However, these sum-integrals all vanish, either because the one-loop 
sum-integral is Z = or X = 0, or because the two-loop sum-integral vanishes: 

1 



fpQ P 2 Q 2 (P + Q) 2 
4 1 



0. 
0. 



T{p} Q P 2 Q 2 (P + Q) 2 
The remaining three-loop sum-integrals have the following forms: 

MiJ ^{PQR}P 2 Q 2 [R 2 Y[(P-Q) 2 y(Q- //)-(//- P)- 

m,j = i 1 



(A.21) 
(A.22) 



pqr P 2 Q 2 [R 2 Y\(P - Q) 2 Y(Q - R) 2 (R - P) 2 

1 



(A.23) 
(A.24) 
(A.25) 



r { P Q} R P 2 Q 2 [R 2 Y[(P - Q) 2 ]i(Q - R) 2 (R - P) 2 
These sum-integrals can be evaluated analytically using methods developed by Arnold and 
Zhai [I3|. The specific integrals that are needed are 
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B Integrals in the Effective Theory 



The effective theory for the scale gT is an Euclidean field theory in 3 space dimensions. Loop 
diagrams involve integrals over 3-momenta. It is convenient to introduce the notation f for 
these integrals. We use dimensional regularization in 3 — 2e dimensions to regularize both 
infrared and ultraviolet divergences. We define the integration measure 

VA 2 V r d 3 ~ 2 * P 



L-fflf 



B.l 



4tt / J (27r) 3 - 2e ' 

If renormalization is accomplished by the minimal subtraction of poles in e, then \x is the 
renormalization scale in the MS scheme. Below, we collect all the integrals that are needed 
to calculate the contribution to the free energy /from the momentum scale gT to order g 5 . 
The nontrivial one-loop integrals that arise in the calculation have the form 
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The specific one-loop integrals that are needed are 
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In ([B.3Q , Iq denotes the derivative of I n with respect to n evaluated at n — 0. 

Some of the two-loop integrals reduce to products of one-loop integrals. The remaining 
two-loop integrals have the following forms: 
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The specific two-loop integrals that are needed are 
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Some of the three-loop integrals reduce to the product of 3 one-loop integrals or to the 
product a a one-loop integral and a two-loop integral. The remaining three-loop integrals 
have the form 
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These integrals are special cases of more general three-loop integrals defined by Broadhurst 
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( 47r )3/2 

Broadhurst derived recursion equations for the integrals Bm and with general arguments 
which can be used to reduce any of the integrals My, N it j and L it j to the basic integrals 
M 0i o and Af ,o, together with simpler 1-loop and 2-loop integrals. The specific integrals that 
are needed in our calculation are 
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We also require the sum of the integrals M^-i and Li,_i, which is simpler to calculate than 
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the individual integrals: 



Mi_i + Li_! = -M 0)0 + 2/i Jj 

- (^(^f P + 0(e)]. (B.30) 
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Figure Captions 

1. One- loop Feynman diagrams for the gluon self-energy. Curly lines, solid lines, and 
dashed lines represent the propagators of gluons, quarks, and ghosts, respectively 

2. Two-loop Feynman diagrams for the gluon self-energy. The solid blob represents the 
sum of the one- loop gluon self-energy diagrams shown in Fig. 1. 

3. One- loop Feynman diagrams for the free energy. 

4. Two-loop Feynman diagrams for the free energy. 

5. Three- loop Feynman diagrams for the free energy. 

6. One-loop Feynman diagrams for the logarithm of the partition function of EQCD. 

7. Two-loop Feynman diagrams for the logarithm of the partition function of EQCD. 

8. Three-loop Feynman diagrams for the logarithm of the partition function of EQCD. 
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